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ELECTRICAL  CONDUCTIVITY  OF  PARTIALLY  IONIZED  GASES* 

ABSTRACT 

A  simplification  of  the  Chapman-Enskog  method  for  the 
calculation  of  the  electrical  conductivity  of  a  multi-component 
partially  ionized  gas  in  a  magnetic  field  is  presented.  The 
calculation  requires  the  inversion  of  a  matrix  which  is  of  the 
order  of  the  approximation,  and  is  independent  of  the  number 
of  species.  The  third  approximation  to  the  electrical  con¬ 
ductivity  is  examined  for  an  electron-ion-neutral  plasma  and 
the  results  are  compared  with  those  obtained  from  the  mixture 
rules  of  Lin,  Resler,  and  Kantrowitz,  and  of  Frost.  It  is 
shown  that  within  the  uncertainties  in  the  experimental 
electron-neutral  cross-section  values.  Frost's  formula  offers 
a  satisfactory  method  of  calculation  for  most  engineering 
applications . 


* 

This  report  has  been  submitted  for  publication  to  the 
Journal  of  the  AIAA. 
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ELECTRICAL  CONDUCTIVITY  OF  PARTIALLY  IONIZED  GASES 
I .  INTRODUCTION 

At  the  present  time,  the  Chapman-Enskog1  method  of  solution 

of  the  Boltzmann  equation  appears  to  offer  the  soundest  theoretical 

basis  for  the  calculation  of  the  electrical  conductivity  of  a 

partially  ionized  gas.  Comparison  of  calculations  made  by  this 

method  with  those  based  on  the  mixture  rules  proposed  either  by 

2  ^ 

Lin,  Resler,  and  Kantrowitz  or  by  Frost  ,  should  provide  a  basis 

for  Judging  the  accuracy  of  these  mixture  rules. 

The  Chapman-Enskog  method  as  extended  by  Hirschfelder, 

4 

Curtiss  and  Bird  for  a  multi-component  mixture  gives  rise  to  a 
set  of  simultaneous  linear  equations  requiring  the  inversion  of 
a  matrix  of  the  order  of  the  number  of  species  multiplied  by 
the  order  o'*  the  approximation.  In  the  limit  of  a  fully  ionized 

C  £ 

gas  it  is  well  known  from  the  work  of  Landshoff  ,  Spltzer  ,  and 

7 

Marshall  that  the  first  and  second  Chapman-Enskog  approximations 
yield  values  for  the  electrical  conductivity  that  are  51  per  cent 
and  98  per  cent,  respectively,  of  the  Spitzer  value.  Therefore 
at  least  the  second  order  approximation  is  required  for  a 
reasonable  level  of  accuracy. 

In  the  Lorentzlan  limit  of  a  weakly  ionized  gas  the  rate  of 
convergence  of  successive  approximations  is  very  sensitive  to  the 
speed  dependence  of  the  electron-neutral  elastic  scattering 
cross-section.  For  the  case  where  the  electron-neutral  collision 
frequency  vgn  depends  on  relative  speed  g  in  accordance  with 
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the  relation  vgn  a  g  ,  the  rate  of  convergence  of  successive 

Chapman-Enskog  approximations  is  illustrated  in  Pig.  1.  (See 

Appendix  B.)  For  an  electron-neutral  collision  cross-section 

tuat  is  linear  with  electron  speed,  the  calculation  must  be 

carried  to  sixth  order  to  obtain  90  per  cent  of  the  final  value. 

For  some  atoms,  the  electron-neutral  elastic  scattering  cross- 

section  is  an  increasing  function  of  relative  speed  in  the 

vicinity  of  thermal  speeds  of  interest,  in  which  case  high  orders 

of  approximation  are  necessary  for  acceptable  accuracy.  The  rate 

of  convergence  for  a  weakly-ionized  gas  calculated  in  the 

Lorentzian  limit  using  experimental  data  for  the  electron-argon 

11 

momentum  transfer  cross-section  is  shown  in  Fig.  2.  For  a 
gas  consisting  of  seeded  combustion  products  where  there  may  be 
of  the  order  of  seven  (or  more)  species  present  in  significant 
amounts,  extensive  calculations  of  the  electrical  conductivity 
to  higher  orders  of  approximation  very  soon  become  Impractical 
for  most  engineering  applications. 

A  simplification  of  the  Chapman-Enskog  method  for  a  multi- 
component  partially  ionized  gas  in  applied  electric  and  magnetic 
fields  is  presented  in  which  the  order  of  the  resulting  matrix 
is  the  same  as  the  order  of  the  approximation  and  is  independent 
of  the  number  of  species.  This  simplification  results  from  being 
able  to  decouple  approximately  the  Boltzmann  equation  for  the 
electron  velocity  distribution  function,  from  the  Boltzmann  equa¬ 
tions  for  the  remaining  heavy  particles  in  the  mixture. 
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The  first  three  approximations  of  the  parallel  electrical 
conductivity  (zero  applied  magnetic  field)  as  a  function  of 
degree  of  ionization  have  been  calculated  for  a  three  component 
gas,  taking  the  electron-neutral  collision  frequency  to  be 
proportional  to  a  power  of  the  relative  electron  speed.*  These 
calculations  are  compared  with  several  proposed  mixture  rules. 

It  is  concluded  that  the  Lin,  Resler,  Kantrowitz  rule  under 
some  circumstances  may  overestimate  the  electrical  conductivity 
by  as  much  as  70  per  cent.  On  the  other  hand,  the  relatively 
simple  mixture  rule  proposed  by  Frost  agrees  with  the  moi  » 
accurate  calculation  to  within  15  per  cent  over  the  complete 
range  of  degree  of  ionization  and  for  a  wide  variety  of  collision 
frequency  models.  Calculations  with  actual  cross-sections 
corresponding  to  cesium-seeded  argon  verify  the  order  of  magnitude 
discrepancy  between  these  two  mixture  rules. 

In  order  to  evaluate  Frost's  proposed  mixture  rule  for 
engineering  applications  we  have  examined  the  spread  in  the 
calculated  electrical  conductivity  resulting  from  uncertainties 
in  the  electron-neutral  collision  cross-sections.  Relying  on 
published  data,  reasonable  upper  and  lower  bounds  for  the 
electron-argon  and  electron-cesium  elastic  scattering  cross- 
section  have  been  constructed  and  these  have  been  used  to  estimate 
the  resulting  uncertainty  in  calculated  electrical  conductivity 
(using  Frost's  formula).  At  temperatures  applicable  to 

*  c 

For  this  special  model,  our  results  agree  with  those  of  Shkarofsky® 
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magnetohydrodynamic  power  generators,  this  spread  in  the  calculated 
conductivity  for  cesium-seeded  argon  exceeds  significantly  any 
difference  between  the  simplified  Chapman-Enskog  third  approximation 
and  the  Frost  formula.  Be  ;ed  on  the  accuracy  of  present  electron- 
neutral  cross-section  measurements,  the  magnitude  of  the  uncer- 

Q  1A 

tainties  we  have  used  appear  to  be  typical  *  .  We  may  therefore 

conclude  that  tnt  mixture  rule  proposed  by  Frost  provides  a 
convenient  and  satisfactory  method  for  calculating  the  electrical 
conductivity  of  partially  ionized  gases. 


II.  SIMPLIFICATION  OF  THE  CHAPMAN- ENSKOG  APPROXIMATIONS  FOR 
CALCULATING  ELECTRICAL  CONDUCTIVITY 

For  purposes  of  brevity,  we  3hall  use  similar  notation  and 

g 

refer  freely  to  the  book  by  Chapman  and  Cowling  .  The  subscripts 
one  and  two  will  be  used  to  denote  electrons  and  ions  respec¬ 
tively,  and  further  subscripts  3,  .  ..,  v  will  denote  neutrals. 

In  accordance  with  the  Chapman-Enskog  approximation  scheme,  one 
seeks  a  solution  of  the  v  Boltzmann  equations  for  the  velocity 
distribution  functions  of  each  species  in  the  form 
fj  =  fj°)[l  +  0 j  +  . ..]  for  j  =  1,  2,  ...  v  where  fj0^  is  a 
Maxwellian  function.  Following  the  development  in  Section  (18.^) 
of  Chapman  and  Cowling,  the  generalization  of  Eq.  (l8.^-10)  for 
the  electron  function  0^(0^)  to  case  of  a  multi-component 
partially  ionized  gas  (as  distinct  from  a  binary  mixture)  is 


where  in  the  present  problem 
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We  have  used  the  synbol  B  for  the  magnetic  induction  (M.K.S. 
units)  and  have  denoted  the  electric  field  in  a  frame  of  refer¬ 
ence  having  the  mean  mass  velocity  of  the  mixture 
cQ  by  E1  =  E  +  cq  x  B  .  The  quantity  I^j(n>g)  denotes  the 
differential  cross-section  for  the  scattering  of  an  electron 
having  speed  g  =  |c1  -  cj  relative  to  a  J-type  particle  into 
the  solid  angle  dn  about  the  unit  vector  ^  .  The  charge 
density  is  denoted  by  pc  and  e-^  =  -e  is  the  charge  on  an 
electron.  The  electron  velocity  is  ,  and 

Ci  =  cx  -  cQ  =  ( 2kT/ ^  . 

Since  Eq.  (l)  Involves  the  unknown  functions  02,  0^,  ...  0y 
in  addition  to  0-,  ,  it  is  necessary  in  general,  to  consider  the 
simultaneous  solution  of  a  set  of  v  Integral  equations  of  which 
Eq.  (l)  is  a  proto-type.  However,  because  the  electron  mass  is 
much  less  than  the  mass  of  any  other  particle  in  the  mixture,  we 
may  take  as  an  approximation  that  the  velocity  of  a  heavy  par¬ 
ticle  is  unaltered  by  electron  collisions,  i.e., 

0j(c')  =  0 j ( cj  for  J  =  2,  3 ,  ...  v  .  With  the  small  electron 


mass  approximation  we  may  also  drop  the  first  term  on  the  right 
hand  side  of  Eq.  (i),and  re-write  the  equation  for  0^  as 


1* 
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By  this  means,  we  are  able  to  de-couple  the  equation  for  the 
electron  velocity  distribution  function  from  the  remaining 
equations.  In  the  small  electron  mass  approximation  the  ion 
current  may  be  neglected  relative  to  the  electron  current  so 
that  the  current  density 


depends  only  on  <t> ^  .  Therefore  to  determine  the  electrical 

conductivity  we  need  concern  ourselves  only  with  the  solution 
of  Eq.  (3). 

In  Appendix  A,  it  is  shown  that  in  the  absence  of  thermal 
diffusion, 

;  (5) 

where  ^  is  a  "generalized"  electric  field,  B  =  B  33  , 

c1  r  cM 

and  C  =  C  -  C  .  In  the  £th  Chapman-Enskog  approximation,  the 
transverse  electrical  conductivities  for  an  electrically  neutral 
plasma  have  the  values 


CT>icr  =  filet 

Wl. 


) 


<PQI 
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where  (£)  is  determined  as  a  solution  of  the  £ 

equations 


N  ^  ~  %  ^rftO  • 


linear 


The  coefficients  Qmm,  are  expressed  in  Appendix  A  in  terms 
of  certain  "bracket  expressions"  which  depend  on  collision 
integrals  ,  defined  by  Eq.  (A23).  These  collision  inte¬ 

grals,  in  turn,  are  determined  by  the  elastic  differential 
scattering  cross-sections  between  electrons  and  other  species 

in  the  mixture.  The  parallel  conductivity  a  =  11m  a-,  ,  and 

B-0  1 

is  therefore  the  same  as  the  electrical  conductivity  in  the 
absence  of  a  magnetic  field. 
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III.  THE  THIRD  CHAPMAN -ENSKOG  APPROXIMATION  (B  =  0) 


Beginning  with  this  section,  we  shall  restrict  our  dis¬ 
cussion  to  the  first  three  Chapman -Enskog  approximations  and 
to  the  case  of  no  magnetic  field;  (the  case  of  B  ^  0  will 
be  discussed  in  a  forthcoming  publication).  For  the  charged 
particle  interactions,  we  shall  employ  the  Rutherford  differ¬ 
ential  scattering  cross-section 


with  the  usual  Debye  length  cut-off.  The  validity  of  this 
approach  is  examined  in  Appendix  C.  Using  Eq.  (8),  the  collision 
integral  for  the  charged  particle  interactions  satisfy  the 
relations 


IS 


IS 


(3) 


Using  the  values  of  the  bracket  expressions  in  terms  of 
the  collision  integrals  provided  in  Refs.  1  and  b,  the  third 
Chapman-Enskog  approximation  for  the  electrical  conductivity 
may  be  written  as 


CttV  a  rfe.1 

a  mt 


* 
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8n,C&+f  ‘f  •  ^  ^  f 

(8/2+26^^+^  m<r"k)X\St?f  (*0 

<6H+*kY$tff-  (JffftlpK  c5»<(* 

and  where  A^  denotes  the  determinant  formed  from  A  by  ^ 

deleting  the  first  row  and  column.  The  qmm  quantities 
are  defined  in  terms  of  the  electron-neutral  collision  inte¬ 
grals  as  follows: 


Equations  (9),  (10),  (ll)  and  (12)  determine  the  parallel  elec¬ 
trical  conductivity  for  arbitrary  electron-neutral  cross-sections. 
To  obtain  the  second  Chapman-Enskog  approximation,  one  deletes 
the  last  row  and  column  from  A  ,  and  similarly  for  the  first 
approximation. 


where 

As 
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IV.  RESULTS  FOR  A  THREE  COMPONENT  PLASMA 

For  the  sake  of  simplicity,  we  shall  consider  a  three 
component  plasma,  and  represent  the  effects  of  electron-neutral 
collisions  by  the  model  Isotropic  differential  cross-section 


2 yn-\ 


where  m  Is  a  parameter. 

Using  Eq.  (13)  one  obtains  the  relation 


and  therefore 


<f*en3n"  , 

^  1 

^  2n3(wf+2w>>)'iwi’'+uwi  +  .  J 

To  examine  the  dependence  of  a(3)  on  degree  of  ioniza¬ 
tion  (n-^/n^)  ,  it  Is  convenient  and  Instructive  to  normalize 
a(3)  with  respect  to  a&dd  ,  where 


(is) 


1- 

<Zi  vti 


(16) 


Equation  (l6)  is  a  mixture  rule  for  calculating  the  electrical 
conductivity  of  a  partially  ionized  gas  proposed  in  Ref.  (2), 


in  which  co 

*  J  ® 


is  the  conductivity  for  a  Lorentzian  gas1  and 

n-  ,  ( 0.582 ~)64(;TU<t£o  (Vlf*  (lB) 

ei  KY^e^A 

g 

is  the  Spitzer-Hclrni  conductivity  for  a  fully-ionized  gas. 

Here,  k  is  Boltzmann's  constant,  T  is  the  absolute  temperature, 
ven(g)  is  the  electron-neutral  collision  frequency,  e  the 
electronic  charge,  €Q  the  dielectric  constant  of  free  space, 
and  A  is  the  ratio  of  the  Debye  length  to  the  impact  parameter 
for  a  90°  collision.  The  physical  basis  underlying  Eq.  (l6)  is 
that  the  resistivity  is  proportional  to  electron  collision  fre¬ 
quency,  and  therefore  to  the  sum  of  electron-neutral  and  of 
electron-ion  collision  frequencies.  Eq.  (l6)  has  the  property 
of  yielding  the  correct  result  in  the  limiting  cases  of  either 
zero  or  infinite  degrees  of  ionization. 

For  a  three  component  plasma,  ven(®)  =  S  >  and 

using  Eq.  (13)>  we  obtain 
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A3  a  convenient  measure  of  the  degree  of  ionization,  we  may 
introduce  the  dimensionless  variable 


where 


) 


yi  _  r(w,|,^/t)rcw- */i\ 
Vs  6(0.582.) 


(20) 

(21) 


Employing  Eqs.  (14),  (15),  (l6),  (l8),  (19)  and  (20),  we  may 


write  Eq.  (10)  in  the  form 

^2.072X 


q-(^_  5TUHX) 
CUi  52(0.5V) 


L+j.smx 

■m ,  ,  • 


ie»+i.%isx 


k„+5.S82X 
X  +l0.78iX 


«<w+X 

a+?oi2x 

W1.S82X 


g,+i.S6ZX 


; 


where 

0(w="  } 

|(2wi+2w*5')')f,.  > 

^=^(wf+2wi+^wl'rttwi+1V2)^  , 
e„-  U^)t . 

4 
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Figures  3  to  6  show  the  dependence  on  degree  of  Ionization 
of  the  first  3  approximations  to  the  electrical  conductivity, 
given  by  Eq.  (20),  for  several  values  of  the  parameter  m 
These  curves  indicate  the  rate  of  convergence  of  successive 
approximations.  For  larger  values  of  the  degree  of  ioniza  'n, 
charged  particle  interactions  dominate.  The  conductivity  is 
insensitive  to  the  parameter  m  ana  convergence  is  rapid.  The 
behavior  at  lower  levels  of  ionization  is  different  however. 

For  m  =  -1  and  for  m  =  0  (Maxwellian  molecules)  convergence 
is  excellent  and  the  third  approximation  seems  to  be  sufficient 
for  all  levels  of  ionization.  For  m  =  1/2  (hard  spheres)  and 
m  =  1  ,  convergence  at  low  ionization  is  less  r^pld.  While  the 
third  approximation  for  m  =  1  is  satisfactory  at  high  levels 
of  ionization,  it  is  not  very  good  at  low  ionization  levels. 

This  behavior,  of  course,  is  merely  a  reflection  of  the  results 
illustrated  for  a  Lorentzian  gas  in  Fig.  1.  Finally,  it  appears 
that  the  mixture  rule  proposed  by  Lin,  Resler,  and  Kantrowitz  may 
in  some  circumstances  over-estimate  the  electrical  conductivity 
by  about  70  per  cent. 
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V.  EXAMINATION  OP  FROST’S  MIXTURE  RULE 

In  1961,  Frost^  proposed  calculating  the  electrical  con- 
ductiA'lty  of  a  partially  Ionized  gas  using  Eq.  (17),  valid 
for  a  Lorentzian  gas,  but  replacing  the  electron-neutral  colli¬ 
sion  frequency  with  the  sum  of  the  electron-neutral  and  electron- 
ion  collision  frequencies.  To  take  into  account  electron- 
electron  interactions.  Frost  introduced  a  modified  expression 
for  the  electron-ion  collision  frequency  having  the  form 

C24) 

in  which 

u  .  AIL  2Tini ei&A  1 2U  \*  I  \'k 

A1  («.«.>»  *  W )  vfT'  > 

is  determined  so  as  to  yield  the  Spitzer-Hclrm  value  of  the 
electrical  conductivity  in  the  fully-ionized  limit. 

In  order  to  examine  the  validity  of  Frost's  procedure,  we 
have  included  in  Figs.  3  to  6  the  results  obtained  using  Frost's 
formula.  At  very  low  values  of  ionization  Frost's  results  seem 
to  be  more  satisfactory,  as  expected,  since  his  formulation  is 
rigorous  in  the  limit  of  zero  electron  concentration.  For  high 
degrees  of  ionization.  Frost's  results  agree  well  with  the 
third  approximation,  and  in  the  fully-ionized  limit  his  results, 
as  they  must  be,  are  slightly  better  than  the  third  approxima¬ 
tion  (by  le3S  than  2%).  Over  the  complete  range  of  ionization 


i 
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and  for  all  the  cross-section  models  examined.  Frost's  results 

and  the  third  approximation  agree  to  better  than  15  per  cent. 

Considering  the  uncertainties  In  actual  cross-seccion  data  (see 

Section  V)  this  agreement  appears  quite  satisfactory. 

In  the  range  of  degree  of  ionization  corresponding  to 
aen 

.01  <  <1  it  is  not  clear  to  what  extent  the  difference 

aei 

between  Frost  and  the  third  Chapman-Enskog  approximation  is 
attributable  to  Frost's  representation  of  the  charged  particle 
interactions,  or  to  a  slow  rate  of  convergence.  To  illustrate 
this  point  we  have  plotted  in  Figs.  7,  8  and  9  for  three  values 
of  m  ,  conductivities  using  Frost's  formula,  but  with 

Hi<Y>,Ki2rS  ; 

and  with 

We  have  also  replotted  the  curves  obtained  with  Eq.  (24)  and  the 
third  Chapman-Enskog  approximation.  The  constants  and  K2 

have  been  adjusted  in  the  same  manner  as  K  .  An  additional 

Q 

comparison  is  provided  by  using  Shkarofsky's  fourth  approxi¬ 
mation  which  is  based  on  a  solution  of  the  Fokker-Planck  equa¬ 
tion  resulting  in  the  same  matrix  as  derived  from  Eq.  (7). 

Shkarofsky's  results  are  not  as  readily  applicable  for  use  with 
experimental  e-n  cross-sections,  as  are  the  Chapman-Enskog 
results . 


The  velocity  dependence  of  Eq.  (25)  follows  directly  from 
the  Coulomb  cross-section  of  Eq.  (8).  For  very  low  degrees  of 
ionization  where  the  electron  distribution  function  i~  unaffected 
by  electron-electron  encounters,  we  would  expect  Frost's  formula 
to  apply  using  Eq.  (25).  The  collision  frequency  of  Eq.  (26) 
has  been  suggested  by  Sodha  and  Varshni11,  to  which  Frost  has 
given  a  single  power  fit,  and  is  based  on  the  Spitzer-HSrm 
numerical  solution  for  a  fully-ionized  plasma.  One  would  expect 
the  use  of  Eq.  (26)  with  Frost's  formula  to  supply  a  close 
approximation  at  high  degrees  of  ionization. 

Figures  7,  8  and  9  show  that  Frost's  formula  with  Eq.  (24) 

lies  between  the  curves  obtained  with  Eqs .  (25)  and  (26)  for  all 

degrees  of  ionization.  This  results  suggests  that  the  velocity 

-2 

power  dependence  g  may  be  viewed  as  a  compromise  choice  which 
achieves  a  reasonable  approximation  to  the  conductivity  for  both 
low  and  high  degrees  of  ionization.  The  Frost  mixture  rule 
appears  to  over-estimate  the  conductivity  at  low  degrees  of 
ionization,  and  to  under-estimate  slightly  at  high  degress  of 
ionization.  The  conclusion  for  low  degrees  of  ionization  is 
supported  by  additional  calculations  in  which  Eq.  (17)  is  used 
with  Eq.  (25),  but  with  K-^  determined  by  the  Rutherford 
collision  cross-section,  rather  than  by  adjustment  to  agree  with 
the  Spitzer-Harm  conductivity. 
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VI.  RESULTS  WITH  EXPERIMENTAL  CROSS-SECTIONS 

Figure  10  shows  Frost's  normalized  conductivity  of  cesium 

seeded  argon  as  a  function  of  temperature,  in  the  range  of 

interest  for  magnetogasdynamic  generators.  Curve  I  makes  use 

12 

of  the  e-Ar.  cross-section  given  by  Frost  and  Phelps  and 

18 

curve  II  is  based  on  Brode's  J  cross-section  for  e-Ar. 

Both  curves  use  Brode's^cross-sectlon  for  e-Cs  .  The  same 
dipping  in  the  actual  conductivity  curve  is  observed,  substan¬ 
tiating  some  of  the  conclusions  drawn  from  the  use  of  an 
electron-neutral  collision  model. 

We  have  attempted  to  estimate  the  possible  spread  in  con¬ 
ductivity  calculations  resulting  from  the  uncertainty  in  experi¬ 
mental  values  of  the  collision  cross-sections.  Figure  11  shows 
the  momentum  collision  cross-section  for  e-Ar.  obtained  by 

14 

different  experimenters.  Based  on  O'Malley's  extrapolation 
for  low  energy  and  the  various  experimental  results  as  well  as 
their  own  work.  Frost  and  Phelp^2have  suggested  an  effective 
electron-argon  cross-section  which  is  shown  on  Fig.  12.  In 
subsequent  calculations  we  have  used  this  curve  as  a  lower  bound 
a.id  Brode's  data  with  extrapolation  at  low  energy  as  an  upper 
bound . 

Figure  11  shows  the  momentum  collision  cross-section  for 
e-Cs  .  The  data  spread  at  low  energy  is  quite  large.  To  determine 
the  corresponding  uncertainty  in  calculated  conductivity  we  have 
chosen  a  reasonable  (low  energy)  upper  bound  guided  by  recent 
theoretical  calculations  made  by  Stone  and  Reltz1^ 
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/p  1  £ 

[0  to  .75  (ev)A/  ]  which  has  been  merged  with  Brode's  data 
[1.6  to  10  (ev)  '  ]  ,  and  a  lower  bound  as  shown  in  Fig.  12. 

Typical  results  of  the  calculations  for  upper  and  lower 
bound  conductivities  of  cesium  seeded  argon,  using  Frost's 
formula  and  the  third  Chapman-Enskog  approximation  are  shown 
In  Fig.  13.  As  expected  from  the  study  of  the  electron-neutral 
collision  model  at  low  Ionization,  Frost's  conductivity  Is 
higher  than  the  conductivity  calculated  via  the  third  approxi¬ 
mation.  As  the  degree  of  Ionization  increases  (temperature 
increase),  the  difference  becomes  less  pronounced . 

The  maximum  difference  between  the  Frost  and  Chapman-Enskog 
conductivities  is  about  25$  and  occurs  at  2000°K.  Over  the 
entire  temperature  range,  the  difference  between  the  Frost  and 
the  Chapman-Enskog  third  approximation  conductivities  is  cf  the 
same  order  or  less  than  the  difference  in  calculated  conductivity 
resulting  from  uncertainty  in  experimental  values  of  the  cross- 

sections.  For  comparison  some  recent  conductivity  measurements 
22 

made  by  Harris  are  included. 
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VII.  CONCLUSIONS 

Figures  3>  ^ ,  5 >  and  6  show  that  the  conductivities  calculated 
by  Frost's  method  and  by  the  third  approximation  to  the  simplified 
Chapman-Enskog  method  differ  by  less  than  13  per  cent  for  a  wide 
range  of  collision  cross-section  models.  Figure  13  indicates  that 
the  uncertainties  in  calculated  conductivity  resulting  from  experi¬ 
mental  cross-section  uncertainties  exceed^or  are  of  the  same  order 
of  magnitude  as  the  differences  resulting  from  these  two  methods  of 
calculation.  Examination  of  experimental  and  theoretical  cross- 
section  data  available  for  elastic  collision  of  electrons  with  a 

variety  of  neutral  species  indicates  that  the  magnitude  of  the 

9  10 

uncertainties  we  have  used  appears  to  be  typical.  *  We  may 
therefore  conclude  that  Frost's  mixture  rule,  at  present,  provides 
a  satisfactory  approximation  for  calculating  the  electrical  con¬ 
ductivity  in  engineering  applications . 
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APPENDICES 

A.  DERIVATION  OP  SIMPLIFIED  CHAPMAN -ENSKOG  APPROXIMATION 
FORMULAE 

Eecause  Eq.  (3)  is  linear  in  <t>^  ,  the  solution  will  have 
the  structure 

"S *%*•"&* 


where  the  unknown  vector  function  D  is  determined  by  inserting 
Eq.  (Al)  into  Eq.  (3)  and  equating  coefficients  of  similar 


gradients.  The  resulting  Integral  equation  for  D  is: 

i\ i  1  W'i  4 


(n) 


The  quantities  B  and  A  are  determined  by  similar  equations, 
but  we  shall  not  consider  these  since  the  dependence  of 
on  E*  (and  thereby  the  electrical  conductivity)  is  provided 
through  D 

For  the  reasons  discussed  in  Ref.  (l)  and  ( 7 ), 

The  quantities  D1  ,  Dg  »  and  D^  are  scalar  functions  of  the 
magnitudes  C-^  and  B  ,  and  are  determined  by  the  three  simul¬ 
taneous  equations 


ttfcoSr  ^  V«Cl}Iij 


(A4-) 


°-f  ^cck- 

tyc^C,- 

jw  ^  f(>  jo  r  ,i  .  t 

°= -*  11 W MW>  ci+ ^d3cc  )c  - 


&*) 


(“0 


Multiplying  Eq.  (A6)  by  B2  and  adding  to  Eq.  (a4)  yields  the 


equation 


|  ||)  y.  to)  f  *2m  \  I 

\  &^Ci4C£)  l^l(Cl^  °i(  )]  C.'- 

Lw^cylc -y^C)  DjCOlQ^-  iQ  Ja  ^ 


for  (D1  +  B2D3)  .  Multiplying  Eq.  (A5)  by  IB  and  adding  to 
Eq.  (a4)  yields  a  single  complex  equation 

V  '^Wi  fi^^LYcc^cVS^cc'K’- 


rCcXcZ^CcK^l^k  ^ 


<**) 
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for  £  =  D.^  +  IB  D2  .  The  electron  cyclotron  frequency  Is 
denoted  by 

•  (AV) 

According  to  Eq.  (18.41-1)  of  Ref.  (l),  (7)  and  Eq.  (4), 
the  contribution  of  D  to  the  current  density  is 


(AlO) 


If  d,  is  decomposed  into  components  parallel  and  perpendicular 

\\  1 

to  the  magnetic  field,  6^  =  d^  +  d^  ,  then 


\(c\  1  Asi 

o-£'W1£+a;|'i&  , 


where 


V 


£ „F  - 

-  ~ 


(All) 


is  a  "generalized"  electric  field,  and  where  B  =  B  B  .  The 
electrical  conductivity  parallel  to  the  magnetic  field  a  ,  is 
determined  by  (D-^  +  B  D^)  and  from  Eq.  (A7),  o  is  inde¬ 
pendent  of  B  .  The  transverse  electrical  conductivities  are 
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determined  in  terms  of 


(Aft) 


Since  reduces  to  a  when  B  =  0  ,  we  need  be  concerned 

only  with  the  solution  to  Eq.  (A8).  Let  us  look  for  a  solution 


of  the  form  of  a  sum  of  Sonlne  polynomials, 

res  Sea , 

rrro  ' 


(AH) 


where  £  is  an  integer  that  defines  the  order  of  approximation. 

In  terms  of  the  coefficients  d^  ,  the  transverse  conductivities 


are  given  by  Eq.  (A13)  as 


JL  X 


-  J) 1- 

»i  V  f €i 


(Aft) 


where  we  have  used  the  orthogonality  property 


(A1 C) 
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Substituting  (Al4)  into  Eq.  (A8),  we  obtain 


,  i  »  •.  (M7 ) 

^c's^trVs&O-^cqsMjyaiA . 


If  we  now  take  the  scalar  product  of  this  equation  with 
S^/2  (^)^i  and  Integrate  with  respect  to  dc^  ,  we  obtain 
the  set  of  4  equations 


VfI=o 


where 


2-  («>+Vi)l 


iwt\,  -fe-  v-" :  'ii:  1 
(it)'*-  w\ 


(.Ml) 


and  where 

[Es&\’T»j|^t1^<oFl{s,vt6lCio+  - 


(M») 
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JW  JM 


•  CKU) 


These  definitions  correspond  to  the  definitions  of  Section  (4.4) 
in  Ref.  (l).  Values  for  the  bracket  expressions  in  terms  of  the 
integrals 

oo 

i\XM  -K1  (S)  /H*\vz 

^  “tf)  ie  *  i  **•©  a- ,  <■* 

li  o  J 


{Kl\) 


where 


Q$V  «  Jo-  An.  , 


(m) 


are  provided  in  Section  (9.6)  of  Ref.  (l)  and  in  Ref.  (4).  The 
quantity  denotes  the  ^educed  mass;  u.^  =  m^/2  and 

U- 12  =  *  *  *  *  ~  • 
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B.  LORENTZIAN  GAS  RESULTS 


In  thi  .  appendix,  we  primarily  summarize  some  results  for 
a  Lorentzian  gas  as  discussed  by  Chapman  and  Cowling^.  The 
formula  for  the  electrical  conductivity  of  a  Lorentzian  gas, 

Eq.  (17)  follows  from  Chapman  and  Cowling's  Eqs .  (10.5-7), 
(9.33-2),  (18.11-5)  and  is  shown  to  correspond  to  the  Infinite 
Chapman-Enskog  approximation  ^gnC00)*  (This  result  has  also  been 

pi 

derived  by  Allis  using  an  expansion  in  spherical  harmonics.) 
The  first  Chapman-Enskog  approximation  to  the  electrical  con¬ 
ductivity  of  a  Lorentzian  gas  may  be  shown  from  Eqs.  (9.81-1), 
(9.8-8),  (18.11-5)  to  be 


3  r^ej  1 


aen^  =  . 

en  lo  n0mn  fi 


1,1 

3“‘l  “13 


(Bl) 


where  is  defined  by  our  Eq.  (A. 23). 

For  the  special  case  of  an  interparticle  interaction  force 
proportional  to  (interparticle  separatlon)"v  ,  Chapman  and 
Cowling  in  Eq.  (10.53-10)  show  that  the  successive  Chapman- 
Enskog  approximations  may  be  written 


_  /m  2  2  2 

e£-n-T  =  1  +  2 —  +  P  (P  t  I,) —  +  .  . .  to  i  terms,  (B2) 
°en' 1 '  q • 1  q(q  +  1) -21 

where  p  =  (v  -  5)/2(v  -l),q=3-2/(v-l)  .  From  Chapman 
and  Cowling's  Eq.  (10.3-8),  the  parameter  v  may  be  related  to 
our  parameter  m  introduced  in  Eq.  (13)  according  to 
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V-rr  -  2”  •  (S3) 

Combining  Eqs .  (14),  (19)  and  (Bl),  we  obtain  for  this  case 

o  (»)  16c  c 

o—TTT  =  W  ^2  "  +  m)  (BM 

en '  ' 

subject  to  the  restrictions  -  |-  <  m  <  |  .  The  results 
plotted  in  Fig.  1  are  based  on  Eqs.  (B.2)  and  (B.4). 
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2.  EFFECT  OF  SMALL  £nA 

The  kinetic  theory  description  of  collisions  between  charged 

20 

particles  has,  until  recently  ,  required  some  sort  of  ad  hoc 
cut-off  procedure  in  order  to  prevent  the  occurrence  of  divergent 
integrals.  The  source  of  this  difficulty  stems  from  the  long- 
range  nature  of  the  Coulomb  interaction  and  the  treatment  of  all 
collisions  as  two-body  encounters.  In  actuality,  the  interaction 
potential  is  essentially  shielded  by  the  collective  behavior  of 
the  particles  at  sufficiently  large  distances.  A  convergent 
theory  is  obtained  by  introducing  some  device  for  ignoring  colli¬ 
sions  with  large  Impact  parameters  in  excess  of  the  Debye  length. 
The  resulting  theory  is  valid  to  order  (lnA)~  ,  where 
A  =  1.24  x  107  (TVn)1//2  is  the  ratio  of  the  Debye  length  to 
the  impact  parameter  for  a  90°  deflection.  (The  quantities  T 
and  n  denote  the  temperature  in  °K  and  the  number  density 
per  m^  respectively.)  In  particular,  the  widely-quoted  value 
of  Spitzer  and  H&rm  for  the  electrical  conductivity  of  a  fully- 
ionized  plasma  is  correct  to  this  order. 

In  Table  (l),  we  have  calculatea  the  values  of  A  and  in  A 
for  a  range  of  temperatures  and  number  densities  which  encompass 
conditions  expected  in  magnetohydrodynamic  generators.  Typical 
values  of  inA  are  seen  to  range  between  4  and  5  which  implies 
an  uncertainty  of  about  25  per  cent  in  the  Spitzer-H&rm  con¬ 
ductivity. 


TABLE  1 

Values  of  A  and  (lnA) 


ne(n 

.-3) 

Te(°K) 

1018 

1019 

1020 

1021 

1022 

1000 

392 

124 

39.2 

12.4 

3.92 

(5.97) 

(4.82) 

(3.67) 

(2.52) 

(1.37) 

1500 

721 

224 

72.1 

22.4 

7.21 

(6.53) 

(5.-41) 

(4.28) 

(3.11) 

(1.98) 

2000 

1110 

351 

111 

35-1 

11.1 

(7.01) 

(5.86) 

(4.71 ) 

(3.56) 

(2.4o) 

2500 

1550 

491 

155 

49.1 

15.5 

(7.35) 

(6.20) 

(5.04) 

(3.69) 

(2.74) 

3000 

2040 

645 

204 

64.5 

20.4 

(7.62) 

(6.47) 

(5.32) 

(4.17) 

(3.02) 

The  region  to  the  left  of  the  heavy  line  In  this  table  Indicates 
roughly  conditions  for  which  the  uncertainty  In  this  theory  is 
less  than  20  per  cent. 

In  about  I960,  Lenard,  Balescu,  Rostoker  and  Rosenbluth^ 
as  well  as  others,  obtained  a  kinetic  equation  which  automatically 
took  into  account  the  collective  behavior  of  a  plasma,  and 
converged  for  large  impact  parameters.  However,  this  so-called 
"fluid -approximation"  diverged  for  small  impact  parameters,  and 


so  it  was  necessary  to  introduce  an  ad  hoc  small  cut-off  limit. 
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A  more  recent  advance  In  the  kinetic-theory  of  charged 
particles  war  made  in  1963  by  Kihara  and  Aono  .  These  authors 
proposed  a  "unified  theory"  for  combining  both  the  close  and 
distant  encounter  contributions  to  the  kinetic  descriptions  of 
collisions,  in  such  a  way  as  to  yield  a  divergenceless  theory 
without  the  need  of  ad  hoc  assumptions.  The  validity  of  this 
new  theory  is  of  order  A”1  ,  and  thus  provides  ample  accuracy 
for  magnetohydrodynamic  generator  applications.  The  region  to 
the  left  of  the  double-line  in  the  preceding  table  indicates 
conditions  for  which  this  new  theory  has  an  uncertainty  of  10 
per  cent  or  less. 

1  Q 

The  theory  of  Kihara  and  Aono  has  been  applied  by  Ikitawa 
to  the  calculation  of  the  electrical  conductivity  of  a  fully- 
ionized  plasma.  The  method  of  calculation  is  Similar  to  that 
of  Chapman  and  Enskog,  and  involves  an  expansion  in  terms  of  Sonine 
polynomials.  The  results  of  Ikitawa  may  be  incorporated  into  the 
framework  of  the  standard  Chapman-Enskog  theory  with  the 
following  identificiation  of  the  bracket  expressions: 


where 


0  0  O 


/-  0  1  '4 
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In  Fig.  14  we  have  plotted  the  ratio  of  Ikitawa's  second 
approximation  to  the  Spitzer-H&rm  value  of  the  electrical  con¬ 
ductivity,  as  a  function  of  InA  .  For  values  of  InA  =  4  , 

the  more  accurate  value  of  the  electrical  conductivity  is  about 
30  per  cent  higher  than  the  Spitzer-H^rm  value. 


j 
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In  1964,  Rynn  reported  a  measurement  of  the  electrical 
conductivity  in  a  fully-ionized  plasma  which  was  20  per  cent 
higher  than  the  Spitzer-Hirm  theoretical  value.  The  probable 
error  of  Flynn' s  measurement  was  stated  to  be  +  10  per  cent. 

At  the  value  of  InA  =  7  for  which  Flynn's  experiment  was  per¬ 
formed,  Ikitawa's  value  of  conductivity  is  about  10  per  cent 
larger  than  the  Spitzer-HSrm  value,  which  places  Ikitawa's 
theoretical  value  within  the  probable  error  of  the  experiment. 
This  comparison  with  experiment  appears  to  provide  some  support 
for  the  new  theory. 

Under  magnetohydrodynamic  generator  conditions,  the  plasma 
is  only  partially  ionized,  and  so  the  magnitude  of  the  correction 
discussed  above  will  be  proportionately  smaller.  When  the 
charged  particle  encounters  are  equally  important  to  electron- 
neutral  collisions,  we  may  expect  a  correction  of  the  order  of 
15  per  cent. 
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